Natural three-qubit interactions in one-way quantum computing 
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We address the effects of natural three-qubit interactions on the computational power of one- 
way quantum computation (QC). A benefit of using more sophisticated entanglement structures 
is the ability to construct compact and economic simulations of quantum algorithms with limited 
resources. We show that the features of our study are embodied by suitably prepared optical lattices, 
where effective three-spin interactions have been theoretically demonstrated. We use this to provide 
a compact construction for the Toffoli gate. Information flow and two-qubit interactions are also 
outlined, together with a brief analysis of relevant sources of imperfection. 
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Measurement based one-way quantum computation 
(QC) P, is an intriguing alternative to the standard 
quantum circuit model and has recently generated enor- 
mous interest for its application to quantum information 
processing (QIP). From the one-way perspective, the sim- 
ulation of quantum gates is performed by adaptive single- 
qubit measurements on suitably prepared multipartite 
entangled resources known as cluster- states. These con- 
stitute a particular class of the more general graph- 
states whose properties have been extensively stud- 
ied 0. Very recently, the basic features of cluster-state 
based QC and a two-qubit quantum search algorithm 
have been experimentally demonstrated 0, 0] . 

However, the realization of cluster-state based quan- 
tum algorithms is often expensive in terms of qubit re- 
sources, an aspect quite detrimental to the efficiency of 
the quantum computation 7., _BJ- An illuminating ex- 
ample is provided by the n-qubit generalization of the 
simple two-qubit searching for a marked entry algorithm 
realized in 5j. In the standard network model for QC ^3] 
this consists of 0(-\/2") oracle-inversion steps [l^ each 
requiring two n-time controlled- NOT (CNOT) gates. For 
n > 3, these can be made from 4(n — 3) three-qubit 
Toffoh gates (C^NOT) 0, which in a cluster-state based 
implementation require 65 qubits 0. A three-qubit ver- 
sion of the algorithm therefore requires ^ 245 cluster 
qubits 01, a number which makes the protocol suscep- 
tible to even small amounts of noise affecting the cluster 
resource. A way to counteract this difficulty is the use 
of more compact designed cluster configurations which 
simulate three-qubit gates. If the universal three-qubit 
Toffoli gate is realized in a compact way, the number 
of qubits and manipulations needed to perform a given 
task will be dramatically reduced. Unfortunately any at- 
tempt in this direction in a cluster-state based scenario 
seems to be destined to failure. The reason is due to 
the underlying two-qubit structure imposed by the effec- 
tive control-phase (CP) gates used in the construction of 
cluster-states Q: a constraint preventing any compact 
natural three-qubit gate is set. It is therefore interest- 



ing to investigate whether other entanglement structures 
are possible for the multipartite entangled resource, pro- 
viding economical configurations which scale better in 
the presence of noise. In this work we describe one such 
possibility, based on recently demonstrated three-spin in- 
teractions in optical lattices We first introduce the 
entanglement structure of the resource and comment on 
simulations of QIP protocols via measurements. Next 
we provide a physically realizable setup for our proposal, 
in the form of a bowtie shaped optical superlattice. Fi- 
nally, imperfections within the model at the entangle- 
ment stages are briefly addressed. Our proposal allows 
the construction of compact configurations for the simu- 
lation of Toffoli gates in one-way QC and opens up new 
possibilities in the search for conducting economical and 
robust-to-noise quantum algorithms. 

The model- We consider a lattice of qubits with the 
bowtie structure depicted in Fig. ^ (a) , where each qubit 
with logical basis {|0),|1)} is initially prepared in the 
state \+) = (1/V2)(|0) + |1)). For every closed trian- 
gle, an entangling operation is applied between qubits 
i, j and k at the vertices equivalent to a control-control- 
phase gate (C2p), S'^'' = 1^'^''^ - 2|lll)yfc(lll|. For 
convenience we denote this operation by O and use 
ai^i (/ = X, y, z) as the Z-Pauli matrix applied to the i- 
th qubit. The physical mechanism which realizes this 
configuration is addressed later. We first describe a way 
to realize information-flow across the lattice. In order to 
create a path for information to be propagated along via 
measurements, it is necessary to remove the influence of 
particular lattice qubits depending on the protocol be- 
ing performed. It is straightforward to check that due 
to the three-body nature of the entanglement, a mea- 
surement in the single-qubit eigenbasis with outcome 
|0) (|1)) destroys (sustains) entanglement between the 
remaining two qubits. On the other hand, by setting 
the qubit to be removed in |0) or |1) before entangle- 
ment is generated across the entire lattice, a path can 
be formed, the choice being dependent on the required 
shape of the path itself. Setting a qubit to |1) gener- 



2 



ates the well-known Ising-type interaction between the 
other two, while setting it to |0) prevents any interaction 
from being generated. Paths of linear cluster-states may 
then be embedded within the lattice and used to prop- 
agate information using techniques borrowed from the 
cluster-state model 0, as shown in Fig.n(a). This tech- 
nique, which initializes the qubits not involved in a spe- 
cific protocol before introducing the entanglement across 
the lattice, puts the removed qubits in an eigenstate of 
(Tz. This vastly reduces the effects of spreading measure- 
ment or environment-induced noise created at the entan- 
gling stage ■ It is easy to see that two-qubit gates can 
be realized in a similar way to the cluster-state model. 
An example is given in Fig. ^ (a) , where measuring the 
bridging qubit (b.q) in the ay eigenbasis simulates the gate 
U = CNOT(l(g)i?j/^)CNOT 0, with R^J^ a single-qubit 
rotation about z on the Bloch sphere by tt/2. 

In addition to the embedded standard cluster-state 
based manipulation of information, the three-spin en- 
tanglement structure can be exploited to construct com- 
pact three-qubit controlled gates using a small number 
of qubits. One example is shown in Fig. ^ (b), where 
we require the logical qubits to propagate away after the 
interaction via O at the central triangle. It is easy to 
see that an enlargement of the basic three-spin trian- 
gle is necessary. In order to give a better insight into 
this, we have extracted the core entangled resource in- 
volved in simulating a Toffoli gate from Fig. ^ (b) into 
Fig. 13 (a). The enlargement can be achieved by mea- 
suring qubits 7 to 10 in the eigenbasis, a method 
similar to that used in the cluster-state model to re- 
move pairs of adjacent qubits. The byproduct opera- 
tion needed to retrieve the original C^P of the central 
triangle between qubits 4, 5 and 6 after enlargement is 
given by U^^ = a% (g) a% ® aff'^QPf^l O CPf^. Here 
sf is the outcome of the measurement of qubit i in the 
ax eigenbasis with sf = (sf = 1) corresponding to 
l+>.: (h). = (1/V2)[|0) - The CP's in U^, resuh 

from C^P not being in the Clifford group 2]. In order 
to show that the enlarged three-spin triangle can be con- 
catenated with the paths propagating the logical qubits 
toward and away from it, we write the CP operator as 
5-^- = _ 2|ll)ij(ll|, which is applied to qubits i and 
i in a triangle when the third qubit k is in |1). Let 
qubits 1, 2 and 3 in Fig.n(b) encode the states \A), \B) 
and \C) respectively. Then, we entangle and measure the 
qubits as follows: {(»'I^^P^){(»lS'''+'^)\+)r\A)i\B)2\C)s, 
where |+)r = ®\=i\+)i and Pf represents the pro- 
jector for a measurement in the ax,i eigenbasis. This 
procedure realizes a two-site cluster-state based prop- 
agation of the logical qubits 0, with qubits 7 — > 
13 left unaffected. Now, let us consider the trian- 
gle enlargement described previously and a subsequent 
information- flow away from it. The entire process is writ- 

ten as (®l£i/^-)56,13^5,1254.11^6,9,8^54059,10^7.854Jx 




FIG. 1; (a): The lattice structure and propagation of quan- 
tum information. Red (blue) dots represent qubits in |1) (|0)) 
and O denotes the three-spin interaction. The dark arrow 
shows the path of information, which is propagated via <jx 
measurements. A two qubit interaction is also shown, where 
the bridging qubit {b.q) should be measured in the Oy eigenba- 
sis Q. (b): A three-qubit interaction and the paths to bring 
information in/out of the interaction region. 

(0f^i5''^+3)|+)T|A>i|B)2|C)3. As = 

(yi,j,k,ljm,n) using the concatenation rules of prop- 
agation [21, we see that the entire lattice can be en- 
tangled and then measurements performed. To com- 
plete this analysis, we show how an arbitary byprod- 
uct operator changes on propagation through the C^P. 

Let ^j^i{ax,ajaz,dj) denote the byproduct operator for 
any measurement pattern M carried out before the 
gate, where aj is the site-label of the logical qubits, 
with values s^^. and s^^ dependent on the outcomes of 
A4 before aj. Upon propagation through, we obtain 

= n'=i(^,cpSa3) xn-=i(^,<^Si'^Sf°^^'"^), 

where Vj is the operator that exchanges label 1 with j. 
A two-qubit byproduct operator is again produced. For 
both C/s^ and C/s^, any CP cannot be propagated triv- 
ially and it is necessary to remove it straight after the 
gate. This can be achieved by applying two-qubit gates 
analogous to the one in Fig. (a) to logical qubits that 
underwent the C^P. In this case, we can measure the b.q's 
in the ay (az) eigenbasis, resulting in a CP (breaking the 
link) between the logical qubits up to local rotations 0, 
thus reducing f/s^^ to local forms again. The simula- 
tion then proceeds as in the cluster-state model until the 
next three-qubit gate occurs. 

Physical Realization- The trapping of alkali atoms such 
as *^Rb in hexagonal two-dimensional optical lattices 
can be achieved using three pairs of counter-propagating 
laser beams (L, £=*=), tuned between the Dl and D2 
line with A = 785 nm and slightly detuned from each 
other. The pairs are in a lin||lin configuration [l^ and 
propagate along y and (y ± \/3x)/2 respectively, pro- 
viding lattice sites with periodicity An appro- 
priate external trapping field is applied in the z direc- 
tion to confine the atoms to the x ^ y plane. Each 
logical qubit can be embodied by the single-atom hy- 
perfine states \a) = |0) = |F = 1,to/ = 1) and 
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FIG. 2: (a): A C^NOT extracted from Fig.[T](b). Qubit |C> 
is the target and O naturally realizes a C^P. (7a;-measurements 
of qubits 3 and 6 realize Hadamard gates before and after O. 
(b): The bowtie lattice structure created by Vojj, where the 
X and y axes are scaled by a = A/2. Dark (light) regions 
correspond to small (large) positive valued potential shifts, 
with the vertices of the superimposed grid representing lattice 
sites. The width wn ~ 2.8 a of a Gaussian beam (used for 
register initialization) is shown. 



\b) = |1) = \F = 2,m/ = 2), with F and m/ the total 
angular momentum of the atom and its projection along 
z respecti vely . They can then be coupled via a Raman 
transition [lj|, using an excited state |e) embodied by an- 
other hyperfine state. We assume the lattice is initially 
loaded with one atom per site, which can be achieved 
by making a Bose-Einstein condensate undergo a super- 
fluid to Mott insulator (MI) phase transition [l^ll5lflH |. 
In a second-quantization picture, the two-species Bose- 
Hubbard Hamiltonian H = -^a=a,bi'^°'^(i,j)^\'^3 ~ 
if/"" Y.i af"? ) + a,^6|a,6, describes the dynam- 

ics of the cold gas of interacting bosons in the peri- 
odic trapping potential considered. Here, (a|) [a — 
a, 6] are the annihilation (creation) operators for atomic 
species a at site i and the summations are taken over 
nearest-neighbor (NN) sites, indicated by (i, j). C/"" 
(jjab _ jjba^ represents the strength of the homo-species 
(hetero-species) on-site repulsive force experienced by the 
atoms. J" is the homo-species tunneling strength be- 
tween sites i and j. For both U and J, effects from atoms 
in NN and next-NN sites respectively are assumed to be 
negligible For a hexagonal periodic potential, we 

can confine the Hamiltonian to a unit-cell of three sites 
in an equilateral triangular configuration [T^ obtaining 
H = i?(o) + V, where = | E»a/3 U°''^ a\l3\(3,ai and 

^ = - E,a + For j; « f/"^ (Vj), 

a,P € {a,b}, V can be treated as a perturbation. In 
these conditions, leaving the subspace M corresponding 
to the prepared MI (with unit filling-fraction) is energet- 
ically unfavorable for the system: its spectrum is gapped 
and states with higher filling-fractions can be adiabati- 
cally eliminated using standard techniques 0]. An ef- 
fective Hamiltonian is obtained with only virtual transi- 
tions to higher population subspaces [13. By switching 
to the pseudo-spin basis |t) = = 1, = 0) and ||) = 



|n" = 0, = 1), this Hamiltonian can be written as [13 
Heff = EU^A,! + Af V| + A^i/^+i + Af ^H^^i + 
A(3)i/^^,^^^, xfE^^^l^^,], where = 

H^j+i = Ei=.,y<^j'^j+i' HL+h3+2 = '^1^1+1^1+2 and 
= Ei=.,,'^j^|+i^j+2- Here A and the A(')'s 
depend on J" and 11°^^ . By varying the laser parame- 
ters. specific parts of H^ff can dominate over the remain- 
der [131 . In particular, we are interested in the terms con- 
taining Aj°', A^^^ and A^^^. Later, we address a suitable 
choice for the physical parameters in H^f / . To realize 
a C^P from Heff, Aq = A3 = — Ai — tt/8 is required, 
where A^ — X^'^^dt. This is possible by tuning J", 
C/"^ and T . However a restriction is imposed by the con- 
dition U / zJ > 5.8, which guarantees the MI regime with 
one atom per site 19j, where z is the number of NN seen 
by a given site. H^ff can be generalized to the hexagonal 
lattice, with each site having 6 NN, so that the model is 
valid for J/U < 0.03. 

In order to produce the bowtie pattern, we can cre- 
ate an optical potential of period > A/2 Here, 
two laser-beams, at angles ±9/2 to a given direction v 
on the X — y plane produce a two-dimensional standing- 
wave in the direction perpendicular to v on the x — y 
plane with period d = A/[2 sin(6'/2)]. Using two pairs 
of lasers, we can produce the periodic pattern Voff = 
Vq — Vi cos(|fc|y) -I- V2 cos(|fc|-\/3x) that offsets the orig- 
inal hexagonal lattice as shown in Fig. |21 (b) . Here Vi 
(V2) is a potential produced by the first (second) pair of 
lasers, Vo = Vx + V2 and = 27r/A. Voff suppresses 
tunneling between specific sites on the lattice according 
to the pattern in Fig. [21(b). During the time evolution, a 
C^P is realized between the sites of closed triangles only. 

The initialization of the register prior to the entangle- 
ment is achieved by applying Raman transitions to all 
lattice sites. These can be activated by standing- waves 
of period a from two pairs of lasers Li and L2, far blue- 
detuned by A from the transition |{a, 6}) |e) and ori- 
entated along the y-axis. All the sites will be located 
at the maximum-intensity peaks and with the atoms 
initially in |a), a rotation of the qubits into the state |-|-) 
can be achieved. In order to perform information-flow 
as described above, individual qubits along the edges 
of a path must be set to |{0, 1}). It is experimen- 
tally feasible to apply a Raman transition to a bunch 
of qubits by addressing them with two lasers of cross- 
section a ~ 10~^^m^. With a Gaussian radial intensity- 
profile, positioning the beams' center between the atoms 
to be addressed, as shown in Fig.|21(b), applies the same 
transition to all the closest surrounding atoms. Thus the 
qubits may be rotated from |+)to |0) or |1) as needed 
using a sort of blurred removal |2ItI | . The state of the cen- 
ter qubit is irrelevant, as it will be disentangled from the 
rest of the lattice. This allows the overlapping of blurred 
removals and the creation of entangled-state subspaces 
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FIG. 3; (a): The fidelity F of the entanglement generation 
against the rescaled time r and the coupling 62. (b): Similar 
to (a) but with the replacement £2 ^ ei , to study imperfect- 
ness in the coupling. 

separated from the rest of the lattice. Other techniques, 
such as using a diffraction limited lens system , would 
also be suitable for this task. Such methods can also be 
used in preparing cluster-states in square lattice configu- 
rations, although single-atom addressing would obviously 
allow more compact gate constructions in both models. 

Imperfect operation- The realization of a single lattice- 
wide C^P should be carried out within the system's 
coherence time. Both A^^^ and A^'^-' roughly scale as 
J^/U^ (J" = J", Vi), with variability upper-bounded by 
J/U ~ 0.03. One can set J and U so that \X^°'^^^'> \ - 0.3 
Hz and A^^'"*) = (A^") can be adjusted by an appropriate 
Zeeman term Here we take j'' ~ 0, J°- — 2 kHz and 

jjaa ^ 2U°'I^ = 120 kHz, possible using Feshbach reso- 
nances |23ji , which correspond to interaction-times within 
the coherence time of this far-detuned configuration j23| . 
Couplings one order of magnitude larger are possible, 
thus lowering the operation time, with the requirement 
that J and U increase by an order of magnitude. Small 
deviations from the desired values of J and U imply slight 
variations of A^*' 's with A^^''*-' becoming nonzero and af- 
fecting the system. In general, as A'^' ~ lOA^'*-' for the 
choices above, we neglect its effect. Thus, the replace- 
ments 

XU) A(^) + €j (j = 0, .., 3), in Heff allow for the 
formal study of imperfect entanglement-generation. The 
imperfect Hamiltonian evolves the initial state |V'(0)) = 
X]/3 7 (5=0 '^p-fsiO)\P, 7, S) of qubits within a triangle into 
the state \ip{t)). To determine the quality of the dynam- 
ics, we compare \tp{t)) to = C^P|V'(0)) using the 
fidelity F = \{->p^%p{t))\'^ . Fig.|S|(a) shows the result for 
ai{0) — l/2\/2. Here 62 has been allowed to vary with 
£0,1,3 = 0. For convenience, we have set |Ao,i,3| = (7r/8)T. 
For r = 1, corresponding to an application of C^P if 
£2 = 0, there is a noticeable fidelity decay against €2- 
Moreover, for £2 > 0, F = 1 at r < 1, which may allow 
a compensation for the fidelity decay by using a shorter 
interaction-time. However this procedure is ineffective if 
£2 is unknown. Fig. O (b) shows another example of the 
effect of imperfect couplings, where we have allowed £1 
to vary, keeping £0,2,3 = 0, a choice due to A'"'^' changing 



more rapidly than A^^^. It is clear that the parameters in 
Hef f should be accurately tuned if the correct entangle- 
ment is to be realized. 

Remarks- We have shown that one-way QC can be per- 
formed on a three-body based entangled resource. Us- 
ing concatenation methods we have described an eco- 
nomic TofFoli gate simulation, fundamental to compact 
algorithm realizations and robust-to-noise QIP. We have 
also analyzed in detail the feasibility of an optical lattice- 
based implementation of our model. 
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